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Abstract 

^ We present a unified method for analysing form factors in B — > ttWi and B — > K*^/ 

d decays. The analysis provides consistency checks on the q 2 and 1/M extrapolations 

'"" i" 1 necessary to obtain the physical decay rates. For the first time the q 2 dependence of 

^ the form factors is obtained at the B scale. In the B — > ttWi case, we show that pole 

c*j fits to / + may not be consistent with the q 2 behaviour of /°, leading to a possible 

factor of two uncertainty in the decay rate and hence in the value of \V u b\ 2 deduced 
from it. For B — > K*j, from the combined analysis of form factors T\ and T2, we find 
the hadronisation ratio Rk* of the exclusive B — > K*^ to the inclusive b — > 57 rates is 
of order 35% or 15% for constant and pole-type behaviour of T2 respectively. 



1 Introduction 



In this paper we describe a method of fitting lattice results for the matrix elements of the 
"heavy to light" decays B — > irlui and B — > K*^ to extract form factors. For B — > irli/i, 
we need the form factor as a function of g 2 , where q is the four-momentum transferred to 
the leptons. The B — > ttWi decay rate can then be predicted, allowing information on the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element V u t, to be extracted once an experi- 
mental measurement is made. For B — > K*^ we need a result at the on-shell point, q 2 = 0, 
where q is the photon four-momentum. This decay has recently been measured by CLEO [1], 
and is a place where physics beyond the standard model, in particular supersymmetry, may 
be tested. 

In both cases, we have lattice calculations performed with heavy quarks at around the 
charm mass which need to be extrapolated to the bottom mass [2, 3, 4, 5, 6, 7, 8, 9]. After this 
extrapolation we typically obtain the form factors close to the maximum squared momentum 
transfer, <7^ ax . This leaves a large extrapolation to be performed to reach low values of q 2 
which dominate the phase space integral for B — > ttWi, or q 2 = which is needed to find the 
rate for B — > K*^. 

Our method is to use the lattice data as fully as possible. We will show how heavy quark 
symmetry makes it possible to obtain information about the q 2 dependence of the form 
factors for B decays. The q 2 behaviour can now be tested, independently of the heavy quark 
extrapolation. Furthermore, we will use kinematic constraints to guide our extractions of the 
form factors. Consequently, in our final results, we have better control over the extrapolations 
and can check their self-consistency. 

2 Form Factors 

Here we give the standard expressions defining form factors for the two decays. 
2.1 B ->• ttWi 

For this decay we need the vector current matrix element between B and tt states: 

(m\uiMB{p))={p + k-q f ^^ f + ( q 2 ) + q , m% ~ 2 ml f (Q 2 ). (1) 

where q = p — k. In the rest frame of the decay products, / + and f° correspond to 1~ and 
+ exchanges respectively. At q 2 = we have the constraint that 

; + (g 2 =o) = /v=o), (2) 

since the matrix element in equation (1) is non-singular at this kinematic point. 

For zero-recoil, where the B and n four-velocities are equal, and q 2 = q^ ax = {ms — 
m v ) 2 = 26.4 GeV 2 , only the time component of the matrix element is nonvanishing. At this 
point, the coefficient of / + vanishes and /°(gmax) alone can be determined. 

The decay rate is dominated by / + since the contribution of f° vanishes for massless 
leptons. 
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2.2 B -> if* 7 

The matrix element of interest is 



{K*(k, V )\sa^b R \B(p)) = J2C;T t (q 2 ), (3) 



where q = p — k as above, r\ is the A'* polarisation vector and 



C\ = 2e„ uXpV »p x kt>, (4) 

c l = nA m \ - m K*) -vq(p + k )m ( 5 ) 

C l = V-qU- o f , (p + k)X (6) 

\ ln B " l K* / 

For an on-shell photon with q 2 = 0, T 3 does not contribute to the B — > K*j amplitude and 
Ti and T 2 are related by, 

T 1 (q 2 =0)=iT 2 (q 2 =0). (7) 

Hence, for B — > K*^y, we need to determine Ti and/or T 2 at the on-shell point. 

The zero-recoil point occurs when q 2 = g^ ax = — mx*) 2 = 19.2 GeV 2 . At this point 
the contributions from T\ and T 3 vanish, so T 2 (g^ ax ) alone can be determined. The expected 
t-channel exchange particle for the T 2 form factor is the 1 + B s i state. 



3 Heavy Quark Symmetry 

We now turn to the predictions of heavy quark effective theory (HQET) for these matrix 
elements in the limit of infinite 6-quark mass [10]. It is convenient to use the four velocities, 
v and v 1 of the mesons, defined by 

p = Mv and k = mv' (8) 

where M is the mass of the initial pseudoscalar (B) and m is the mass of the final vector 
(K*) or pseudoscalar (tt) meson. We also use the variable 

M 2 + m 2 - q 2 

LO = v ■ v' = 9 

2Mm v ; 

so that zero-recoil occurs at lo = 1. 

The zeroth order HQET predictions for the B — > nWi form factors are: 

n<f,z) = f (Vm-^), Co) 

fx 



f°(q 2 ,x) = ^((l-^fH-^-^)). (11) 

where x = m/M. The functions ^f 2 depend on the light degrees of freedom, but do not 
depend on the heavy quark mass. They are also independent of the gamma-matrix structure 
of the current. 
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For B — > K*~/ the corresponding form factor predictions are [11]: 



Tl ( q \x) = ^(^H-^), (12) 



iT 2 {q\x) = ^((l-^H-^-^H), (13) 

The expressions are identical in form to those for B — > irli/i, but the functions 6\ 2 are 
different since they depend on different light degrees of freedom. We note in passing that 
the same two Q\ 2 functions govern the HQET prediction for the V and A\ form factors in 
heavy pseudoscalar to light vector semileptonic decays proceeding via the left-handed vector 
current. 

These predictions will be subject to perturbatively calculable renormalisation by strong 
interactions and corrections from terms with higher powers of 1/M. However, they provide 
us with a starting point for our fitting procedure. 

For a fixed value of ui, the predictions in equations (10), (11), (12) and (13), give simple 
scaling laws in the heavy quark limit, M — >■ oo (x — >■ 0). In particular, at q^ ax we obtain 

/ + (<?Lx) ~ , T^gLx) ~ 

The point with q 2 = does not correspond to fixed lu as M — >■ oo. This means we have to 
guess the q 2 dependence of the form factors if we are to extract scaling behaviour at q 2 = 0. 
Pole dominance ideas suggest that, 

^ - d-,vV (15) 

for / = /+, /°, 7\, T 2 , where Mf is a mass that is equal to M plus 1/M corrections and rif is 
a power. Since 1 — q^^/M 2 ~ 1/M for large M, the combination of heavy quark symmetry 
and the form factor relations at q 2 = implies that Uf+ = np + 1 and = riT 2 + 1. So, if 
we fit / + or T\ to single pole forms, then f° or T 2 should be constant in q 2 . To have a single 
pole form for f° or T 2 necessitates double pole (or "dipole") forms for / + or T\. These two 
types of behaviour correspond to 

t+(a\ t ( c\\ ) M~ l l 2 single pole 

/ (0) or TWO) ~ < , , o/ 9 n n , , • (16) 
w v ; [ M~ 3 / 2 double pole v ; 

In passing from pole/constant to dipole/pole behaviour we are allowing f° for B — > ttWi 
or T 2 for B — > K*j to acquire curvature in q 2 . As we will see, both f° and T 2 have weak 
dependence on q 2 for the measured points: in this case constrained multipole fits, for example 
tripole/dipole, to these form factors are nearly indistinguishable in the region of q 2 between 
and (/max, since one can compensate higher powers by changes in the fitted mass parameter. 
For this reason, we will show results for pole/constant and dipole/pole fits only. 



4 Fitting Lattice Data 

From lattice calculations with propagating (rather than static) quarks, we can obtain matrix 
elements for heavy quarks around the charm mass over a range of q 2 straddling q 2 = 0. In 
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Figure 1: Labelling of quark hopping parameters for three point correlator calculation. 

extracting form factors from these matrix elements, we can reach g^ ax for f° and T 2 only. 
We need to extrapolate in the heavy mass M to scales around the b quark mass. Such scaling 
in M is simple for fixed u, but, at the B scale, produces a range of q 2 values near g^ ax and 
far from q 2 = 0. We are therefore faced with a large extrapolation to q 2 = 0. 

The results described below come from 60 SU(3) gauge configurations generated by the 
UKQCD collaboration on a 24 3 x 48 lattice at (3 = 6.2 in the quenched approximation. The 
O(a) improved Sheikholeslami-Wohlert [12] action was used for fermions, with "rotated" 
fermion fields appearing in all operators used for correlation function calculations [13]. The 
inverse lattice spacing determined from the string tension is a -1 = 2.73(5) GeV [14]. 

Three point correlators of the heavy-to-light two fermion operators with a heavy pseu- 
doscalar meson (the "5" meson) and a light pseudoscalar or vector meson were calcu- 
lated, as illustrated in figure 1. Matrix elements were extracted from these correlators 
by the method detailed in [9, 15, 16, 17]. Four heavy quark hopping parameters, = 
0.121,0.125,0.129,0.133, were used. For the propagator connecting the current opera- 
tor to the light meson operator, two kappa values, k s = 0.14144,0.14226, were available. 
The subscript s is for strange: these kappa values straddle that for the strange quark, 
0.1419(1) [18]. For = 0.121,0.129, we used three light "spectator" hopping parameters, 
m = 0.14144,0.14226,0.14262, and for n h = 0.125,0.133 we used m = 0.14144 only. The 
critical hopping parameter at this f3 is n cv \i = 0.14315(1) [18]. 

The lattice calculations were performed with the heavy meson spatial momentum of 
magnitude or 1, in lattice units of 7r/12a. The momentum injected at the operator insertion 
was varied to allow the modulus of the light meson spatial momentum to take values up to 
^/S in lattice units (although some of the momentum choices were too noisy to be used in 
fits). We refer to each combination of light and heavy meson three-momenta as a channel 
with the notation |p| — > |k| in lattice units (for example — >■ 1 or 1 — > 1± where the subscript 
_L indicates that p and k are perpendicular). 

The results below have been obtained using uncorrelated fits for the extrapolations in 
the heavy quark mass and in q 2 . The extraction of the form factors from the three-point 
correlation function data used correlated fits [9, 15]. Statistical errors are 68% confidence 
limits obtained from 1000 bootstrap samples. 

To make the best use of HQET, we pick momentum combinations which keep lu constant 
or nearly constant as the heavy mass varies. This will allow us to scale linearly or quadrati- 
cally in 1/M from the charm to the bottom mass scale and gives us the form factors for the 
B decays as a function of to. Now we convert to to q 2 and fit to assumed forms for the q 2 
dependence at the B scale, consistent with the relations / + (0) = /°(0) and Ti(0) = zT 2 (0). 
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momenta 








P 




k 




0.121 


1 25 


1 29 


1 33 


^ave 










1 nn 
1 .uu 


1 nn 
1 .UU 


1 nn 
1 .UU 


1 nn 
1 .UU 


1 nn 
1 .UU 







1 




L.oo 


L.oo 


1 'X'X 
L.oo 


L.oo 


1 .oo 





V2 


i fin 
l.oU 


1 fin 
l.oU 


i fin 
l.oU 


1 fin 
l.oU 


i fin 

1 .UU 




1 







1 04 

J. . V7t: 


1 05 




1 09 


1 nfi 

1 .UU 




1 




u 




1.38 


1.40 


1.42 


1.45 


1.41 




1 


V2 ± 


1.66 


1.68 


1.70 


1.74 


1.69 



Table 1: Values of lu and their averages for B — > ttWi with various values and k s = Ki = 
0.14144. We average over values of the light meson momentum k when possible, but for the 
heavy meson only two momenta are available (p = (0,0,0) and (1,0,0) in lattice units of 
ir /12a). 



momenta 








P 




k 




0.121 


0.125 


0.129 


0.133 


^ave 








1.00 


1.00 


1.00 


1.00 


1.00 







1 


1.22 


1.22 


1.22 


1.22 


1.22 




1 





1.04 


1.05 


1.07 


1.09 


1.06 




L 




U 


1.27 


1.28 


1.30 


1.33 


1.30 



Table 2: Values of lu for B — ^ K*j for various momentum channels and values, with their 
average for each channel. Calculated with m = nix* for k s = 0.1419 and Ki = 0.14144. 



Since 



LU = V ■ V 



p k 



Mm 



(17) 



we can select channels where p • k = 0, so that 



_ E m , p^ 



(18) 



We see that lu is independent of M when the heavy meson is at rest and the light meson 
momentum is fixed. When |p| = 1 in lattice units, the change in lu is only about 6% in a 
given channel for our range of heavy kappa values. Given that this error is comparable with 
others from discretisation effects, quenching and so on, this should not preclude inclusion of 
these channels in fitting. We take the actual lu as an average of the u/s for the four heavy 
kappas. In tables 1 and 2 we show the channels used and the corresponding lu values. 



4.1 B — > ttIvi Results 

For this exploratory study, not enough spectator quark kappa values were available for 
all heavy quark kappa values to allow reliable chiral extrapolations. Therefore, our aim in 
this section is not to quote realistic values for the form factors governing the B — > ttWi decay, 
but to illustrate the main features of the analysis method proposed in this paper. All results 
presented in this section are for the 0.14144. 
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Figure 2: Linear and quadratic scaling of 
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0.2 



0.8 
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and f° for the channel 0—^1. The perturbative 
value, 0.83 [19], was used for the vector current renormalisation constant, Zy. M is given in 
units of GeV. The dotted lines show the statistical error on the linear fits. 



We first scale the form factors measured at fixed u> to the B mass, using both linear and 
quadratic dependence on 1/M to extrapolate: 

/+e/v/M = M 1 + r) \ line " («) 

l 7+ (1 + it + J quadratic 

/ — f 7o f 1 + %) linear 
/WM = I ) t f ) , (20) 
[ 7o (l + f + ^) quadratic 

where 9 comes from the leading logarithmic factors and is chosen to be 1 at the B mass, 

with /So = 11 in the quenched approximation and Aqcd = 200 MeV. An example of the 
scaling is shown in figure 2 for the channel 0—^1. We quote results from the linear ex- 
trapolation, but, as an indication of the systematic error, incorporate the difference between 
linear and quadratic extrapolations in quadrature with the statistical error. We have not 
incorporated any systematic error due to ambiguities in the determination of a -1 . We note 
in passing that taking the chiral limit would also increase our errors. 

In figure 2 we show also the form factors in the static limit obtained from [20]. The dis- 
crepancy between the extrapolation and the static point could be due either to discretisation 
errors for the propagating quarks or to excited state contributions in the static results. The 
limited agreement gives some confidence in our extrapolations. 

The extrapolated form factors were fitted separately to different q 2 dependences: / + to a 
dipole and a pole, f° to a pole and a constant. Momentum channels 0—^1,0—^ a/2, 1 — > 0, 
1 — > lj_, 1 — >■ V2 ± and, for f° only, 0—^0, were used. The fits are illustrated in figure 3. 
The dipole/pole combination for / + , f° is favoured over the pole/constant or indeed over 
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Figure 3: Unconstrained dipole/pole fits (solid lines) and pole/constant fits (dashed lines). 
The dotted lines mark q 2 = and The x 2 /dof for the fits were comparable save for 

the constant fit to f° which was about 4 times greater. 



fitting both form factors to single poles, because this combination comes closest to obeying 
the relation /+(0) = /°(0). 

Dipole/pole and pole/constant fits constrained to satisfy / + (0) = /°(0) are shown in 
figure 4 with numerical values given in table 3. Again the dipole/pole fit is preferred over 
the pole/constant fit. 

Constrained fits of f +,0 (q 2 ) were also carried out for each of the four k^s using momentum 
channels 1 — > — 1, 1 — >■ 1 and 1 — > (1,1,0) in addition to those used above. The resultant 
interpolations to f(q 2 = 0) (we drop the superscripts on /+ and f° when they are constrained 
to agree at q 2 = 0) were scaled appropriately to the B mass according to equation (16). These 
are the "burst points" shown in figure 4, with values given in table 3. The agreement between 
the two methods shows the self consistency of using dipole/pole or pole/constant constrained 
fits. 

In previous lattice calculations [2, 3] the / + form factor was determined using the — >■ 1 
momentum channel and assuming a q 2 dependence given by the exchange of the B* resonance 
in the ^-channel. We will refer to this as the "pole procedure" below. Here we use the physical 
B* mass (5.46 GeV) for the pole since this corresponds roughly to the quark kappa values 
employed. This procedure gives / + (0) = 0.43 ± 0.02, in reasonable agreement with the 
constrained pole/constant fit. 

The results of table 3 and figure 4 suggest that a constrained dipole/pole fit is more 
appropriate than a pole/constant fit. The result of using such a fit for /+ is to reduce the 
value for the decay rate from that using a pole behaviour for / + , as shown in table 4. As a 
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1 2 3 1 2 3 



Figure 4: Constrained dipole/pole (left) and pole/constant (right) fits to / + and f° at the B 
mass. The burst point, is from /(0)'s obtained from constrained fits at the four k^s scaled 
appropriately to the B mass. The shaded bands mark the error on /(0) obtained from the 
fits. The dotted line marks g^ax- 



consequence, the extracted value of the quark mixing parameter \V u b\ 2 changes by a factor 
of two. 

4.2 B -> K*j Results 

The lack of spectator quark kappa values is less of a problem here than for B — > irlvi. It has 
been shown that the form factors' dependence on the spectator quark mass is mild [9]. As 
in reference [9] we interpolate k s to the physical strange kappa value, 0.1419(1), but neglect 
any dependence of the form factors on the light spectator quark mass. Therefore we present 
results only for ki = 0.14144. We then scale to the B mass using a similar procedure to that 
explained for B — > irlvi. 

For T\ we use the following momentum channels: — >■ 1, 1 — >■ and 1 — > 1±. For T2 we 
can use the — > channel in addition to these. We do not use the channels — > \/2 and 
1 — > V2 ± as both the statistical and systematic errors on the points are large. We fit 7\ and 



fit type 


burst pt. 


/(0) 


m f+ 


m f0 


X 2 /dof 


dipole/pole 
pole/const 


0.23 ±0.02 
0.45 ±0.02 


o.24i8;8| 

0.49 ±0.02 


5.7±0.1 
5.6 ±0.3 




1.5 
6.1 



Table 3: Constrained dipole/pole vs pole/constant fits at the B mass for B — > ttWi with 
k s = Ki = 0.14144. The burst point, see figure 4, is from /(0)'s obtained from constrained 
fits at the four k^s scaled appropriately to the B mass. Fitted masses are in GeV. 
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fit type 


T(B - 




v ub 


2 10 12 s -l 


dipole/pole 




c Q+1.2 
°-°-1.0 




pole/const 




13.0 ± 1.6 


pole procedure 


10.2 ±3.6 



Table 4: Decay rates calculated using different forms for f + (q 2 ). The result in the last row 
has been obtained using the — >■ 1 momentum channel and a single pole form for / + , as 
explained in the text. 
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it — interpolated then scaled 
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0.0 
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0.6 — 



0.4 



0.2 



0.0 



Tj.Tg, at B mass 



□ ,o - fitted points 

— interpolated then scaled 




Figure 5: Constrained dipole/pole (left) and pole/constant (right) fits for T\ and T 2 at the B 
mass. The burst point, is from T(0)'s obtained from constrained fits at the four /c^'s scaled 
appropriately to the B mass. The shaded bands mark the errors on T(0) obtained from the 
fits. The dotted line marks g^. 



T 2 simultaneously, imposing Ti(0) = iT 2 (0), with pole/constant or dipole/pole forms. These 
fits are shown in figure 5 with numerical values given in table 5. 

Alternatively, for each Kh, we fit T\ and T 2 as functions of q 2 , with Ti(0) = ^(0) imposed, 
adding channels 1 — » — 1 and 1 — > 1 to the above list. This gives T(q 2 =0) (we let T(0) denote 
Ti(0) = iT 2 (0)). These points are extrapolated to the B scale according to equation (16). 
The q 2 = point is the burst point plotted in figure 5 and referred to in table 5. 

As can be seen in figure 5 and table 5, fitting form factors using the two methods de- 
scribed above produces consistent results. The burst points should agree with the results of 
reference [9]. For dipole/pole fits, this is indeed the case, but for pole/constant fits a much 
lower value was found for T(0) in [9]. We believe this difference arises from our imposition 
of the constraint 2^(0) = iT 2 (0), which has reduced the curvature in the values of T(0) as a 
function of 1/M. In [9] the curvature produced a low value of the form factor at q 2 = at 
the B scale, which is inconsistent with the value of T 2 (g^ ax ) if T 2 is constant. 

We have checked that higher power pole fit combinations, for example tripole/dipole, 
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fit type 


burst pt. 


T(0) 


m Tl 


m T2 


X 2 /dof 


R K * 


dipole/pole 
pole/const 


0.14 ±0.03 
0.25 ±0.02 


1 5 +0 - 07 
u - i0 -0.06 

?fi+ 02 
U - ZD -0.01 


r O+0.6 
5-O-0.5 

5.4 ± 0.1 


7+3 
' -2 


0.4 
0.9 


(13^)% 
(35^)% 



Table 5: Dipole/pole vs pole/constant fits for B — > K*j for k s = Ki = 0.14144. The burst 
point, see figure 5, is from T(0)'s obtained from constrained fits at the four k^s scaled 
appropriately to the B mass. Fitted masses are given in GeV. Rk* is the hadronisation 
ratio defined in the text, where a value of mi, = 4.65 GeV [21] is used. 



give results consistent with the dipole/pole fit. 

In contrast to the B — > ttWi case, our final results do not allow us to favour any particular 
fit combination. Hence, from this study, we can conclude only that T(0) lies in the range 0.09 
to 0.28. This situation is common to current lattice calculations of B — > K*j form factors [6, 
7, 9]. With higher statistics and smaller lattice spacings (allowing larger spatial momenta and 
hence smaller q 2 values), the method proposed here should be able to differentiate between 
q 2 behaviours and eliminate the present uncertainty. 

For comparison with experiment it is useful to consider the hadronisation ratio, given up 
to 0{l/m 2 ) corrrections by [21] 



R 



K' 



Y(B K*j) 




\m b J 



m. 



T(0) 



(22) 
(23) 



in which many of the theoretical uncertainties in relating the form factor to the branching 
ratio cancel. 1 We report values for Rk* in table 5. The experimental value is Rk* = 
(19 ± 13)% [1, 23]. It appears that the dipole/pole fit is in better agreement with the 
hadronisation ratio data, but the errors in both the fit and the experimental measurement 
are large. 



5 Conclusions 

We have presented a unified method for analysing form factors in B — > ttWi and B — > K*j 
decays. We have extrapolated as many momentum channels as possible at nearly fixed lu to 
give the maximum information on the q 2 dependence of the form factors at the B scale. For 
the first time, we have been able to fit this q 2 dependence directly. This analysis, combined 
with the procedure of extrapolating to q 2 = and then scaling to the B mass, provides 
consistency checks on the extrapolations necessary to obtain the physical decay rates. The 
imposition of the constraint at q 2 = in all our fits is novel and provides a framework for 
the best utilisation of the lattice data. 

Our main conclusions are that for the B — > ttWi case, pole fits to / + may not be con- 
sistent with the q 2 behaviour of /°, leading to a possible factor of 2 uncertainty in future 
determinations of \V u t,\ 2 from this decay. For B — > K*j, with the current data, no particular 

1 The theoretical prediction for this ratio may be subject to long-distance effects [22]. 
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q 2 behaviour at the B scale is favoured. The resulting uncertainty in T(0) is large, limiting 
the usefulness of measurements of this decay for constraining the standard model or new 
physics. 

The problems associated with the q 2 extrapolation will be generic to the study of any 
decay of a B meson into a light meson plus leptons or a photon. Direct simulations of the B 
meson on the lattice, using NRQCD or the heavy Wilson action (as proposed by the Fermilab 
group), will face the same problems, if they are limited to a small range of q 2 around g^. 

These considerations do not modify previous results for charmed meson semileptonic 
decays, because the q 2 range available from lattice simulations covers the physically accessible 
range. Pole and dipole fits for / + agree over the physical range of q 2 for charm decays. 

The analysis method presented here, together with measurements closer to q 2 = 0, will 
provide an unambiguous determination of the q 2 behaviour of the form factors. In this study, 
we could not use the 0, 1 — >■ (and 0, 1 — >■ \/2 for B — > K*j) channels which were already 
too noisy. Future lattice simulations with higher statistics and smaller lattice spacings at 
fixed physical volume, will allow more momentum channels to be used, thereby increasing 
the usable range of q 2 . 
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